Fluctuations of the vortex line density in turbulent flows of quantum fluids 
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We present an analytical study of fluctuations of the Vortex Line Density (VLD) < 
5C{uj) 5C{—lj) > in turbulent flows of quantum fluids. Two cases are considered. The 
first one is the counterflowing (Vinen) turbulence, where the vortex lines are disordered, and 
the evolution of quantity C(t) obeys the Vinen equation. The second case is the quasi-classic 
turbulence, where vortex lines are believed to form the so called vortex bundles, and their 
dynamics is described by the HVBK equations. The latter case, is of a special interest, 
since a number of recent experiments demonstrate the w -5 / 3 dependence for spectrum VLD, 
instead of uj 1 ^ 3 law, typical for spectrum of vorticity. In steady states the VLD is related 
to the normal velocity as £ — (p^/ps) 2 v 2 for the Vinen case, and £ = |V x v„| /k for 
rotating vortex tubes. In nonstationary situation, in particular, in the fluctuating turbulent 
flow there is a retardation between the instantaneous value of the normal velocity and the 
quantity £. This retardation tends to decrease in the accordance with the inner dynamics, 
which has a relaxation character. In both cases the relaxation dynamics of VLD is related to 
fluctuations of the relative velocity, however if for the Vinen case the rate of temporal change 
for £(t) is directly depends on Sv ns , for the HVBK dynamics it depends on V x 5v ns . As a 
result, for the disordered case the spectrum < 5£(co) 5£(—lu) > coincides with the spectrum 
w -5/3 j n cage Q £ ^ e b unc u e arrangement, the spectrum of the VLD varies (at different 
temperatures) from a; 1 / 3 to a; -5 / 3 dependencies. This conclusion may serve as a basis for 
the experimental determination of what kind of the turbulence is implemented in different 
types of generation. 
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I. INTRODUCTION 



The problem of modeling classical turbulence with a set of chaotic quantized vortices is the 
lottest topic in modern studies of vortex tangles in quantum fluids (see e.g., recent reviews articles 
lj] j Q] > EJ ) • The most common view of quasi-classical turbulence is the model of vortex bundles. The 
point is that the quantized vortices have a fixed core radius, so they do not possess the important 
property of classical turbulence - stretching of tubes - which is responsible for the turbulence energy 
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cascade from the large scales to the small scales. The collections of the near-parallel quantized 
vortices (vortex bundles) do possess this property, so the idea that quasi-classical turbulence in 
quantum fluids is realized via vortex bundles of different sizes and intensities (number of threads 
) seems to be quite natural. 

Recently two numerical evidences of the vortex bundles structures were obtained. Thus, in 
one numerical work, [J], in 2048 3 simulation of quantum turbulence within the Gross-Pitaevskii 
equation, authors observed nonuniform structures. The authors claimed that "the visualization of 
vortices clearly shows the bundle-like structure,which has never been confirmed in GPE simulations 
on smaller grids." In the other numerical work [ij], (J, the authors studied the evolution of the vortex 
structures (at zero temperature) on the basis of the Biot-Savart law. They also observed structures 
reminiscent of the field of vorticity in classical turbulence (see e.g., Q]). 

As for experimental confirmation, there are not so far strong evidences of the bundle structure. 
On the contrary, there are experimental results, which would seem to refute the idea of bundles. 
Thus, in experiments by Roche et al. g], and by Bradley et al.jj]], it was observed that the 
spectrum of the fluctuation of the VLD C is compatible with a —5/3 power law. This contradicts 
the idea of the bundles structure, since the spectrum of the vorticity (and, correspondingly, of the 
VLD C (via Eq. (HI))) should scale as 1/3 power law. An explanation was offered by Roche and 
Barenghi [lo|. The authors considered the VLD C to be decomposed into two components. The 
one,, smaller, part consisting of the polarized component, is responsible for the large scale turbulent 
phenomena, whereas the other, disordered, part evolves as a "passive" scalar, thereby taking the 



—5/3 velocity spectrum. In paper [ll|], the authors performed the direct numerical simulations of 
the "truncated HVBK" model. They confirmed the existence of the u; _5//3 spectrum for fluctuation 
of the VLD, however, for the larger temperature this dependence became more shallow (probably 
reaching W 1 / 3 ), as it should be for the classic turbulence. 

We present an analytical evaluation of the spectrum of fluctuations VLD < 5C{uS) 8C{—uj) >. 
Two cases are considered. The first one is the counterflowing (Vinen) turbulence, where vortex 
lines are disordered and dynamics of quantity C(t) is governed by the Vinen equation. The second 
case is the quasi-classic turbulence, where the vortex lines are believed to form the so called vortex 
bundles, and the dynamics of VLD obeys the HVBK equations. In steady states the VLD is related 
to the normal velocity as £ = (p'y/ps) 2 v\ for the Vinen case and C = |V x v n | /k for the rotating 
vortex tubes (notations are standard, see e.g. fl2]). In nonstationary situation, in particular, in the 
fluctuating turbulent flow there is a retardation between instantaneous value of the normal velocity 
and the quantity C(t). This retardation tends to decrease, according to the inner dynamics, which 
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has a relaxation character. In both cases, the relaxation of 5C(t) is related to fluctuations of 
the normal velocity 5v ns velocity. If, however, for the Vinen disordered turbulence, the rate of 
temporal change for £{t) is directly depends on 5~v ns , for the HVBK dynamics it depends on the 
quantity V x 5v ns . In addition, the relaxation mechanisms, and, consequently, times of relaxation 
are different. The factors, outlined above lead to different formulas for spectra < 5C{uj) 5C(—oj) > 
and their dependence on temperature. 



II. VINEN EQUATION CASE 

Let us study reaction of the vortex line density, in fluctuating flow of normal velocity supposing 
that the dynamics of £{t) obeys the Vinen equation 



^= a v \v ns \C^ - p v £ 2 . (1) 
equation ([!]) was derived phenomenologically jnitially for pure counterflowing superfluid helium 



13Hla] . Attempts to derive it an analytic form 



17K2fJ] demonstrated that this equation is seemingly 



valid for any non-structured turbulence. Under term "non-structured turbulence" we understand 
the vortex tangle, which consists of closed vortex loops of different sizes, uniformly distributed in 
space. It differs, for instanse, from the turbulent fronts in rotating fluids, which deals with the lines 
terminating on lateral walls. It also differs from the the mechanically excited turbulence, which is 
believed to consists of the so called vortex bundles composed of very polarized vortex filaments. 
Beside of the counterflow turbulence, the "non-structured turbulence" is generated by intensive 
sounds (both by the first and second). The case of vortex tangles, which appear also after the 
quench due to the Kibble-Zurek mechanism, or by the proliferation of vortices when approaching 
the critical temperature. 

Our goal now is to study the stochastic properties of C(t), when v n fluctuates with the a 
given spectral density (5v n (u;)5v n (— to)) = J(uj). Further, for simplicity, we will study the pure 
counterflowing case in sense that the net flow is absent, j =p n v n + p s v s = 0. Then, the average 
value £oof the vortex line density is related to a relative velocity v ns = v n — v n by the usual 
relation 



£ n = ^Zv 2 = -y 2 v 2 (2) 

a2 ns I v ns' \^) 

Py 



To take into account the fluctuations, we put 

£ = 7 2 v ns + v n = v «0 + S\ n . (3) 
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From condition j =0 it follows 



Pn> c P c / a \ 

v s = v s0 ov n , v„ s = v„ s0 H ov n . (4) 

Substituting Eqs. ([3]) and (H]) into the Vinen equation ([1]) we arrive at 

= ay —£l /2 5v n - p v ^-£ S£. (5) 
at p s 2 

Equation ([5]) shows, that the evolution of fluctuating part of the vortex line density 5C bears the 

relaxation-type character with a characteristic time Ty = 2/ {(3yCo) and with a "force" proportional 

to 5v n . It allows to express the spectrum of VLD (5£(uj)5£(—uj)) via the spectrum of a normal 

velocity. We accept further the usual in theory of turbulence relationship k = uj/v n between the 

wavenumber k and the frequency uj (v n is the mean flow). In the Fourier transforms Eq. ([5]) takes 

a form 



iuj5£(uj) = ay — C^j 2 5v n {uj) — f3y-£o5£.(uj), (6) 

Ps ^ 

therefore the spectrum (S£(u)S£(— u)) is 

4(a v /(3y) 2 (f ) 2 £ {6v n (u)6v n (-uj)) 
(5£(u)5£(-u)) = i + (a,Tv)2 ^ 

Relation ([7]) shows that for small frequencies, uj < 1/r, the spectrum of the VLD reproduces 
the spectrum of fluctuations of the normal component, and if the Kolmogorov-type turbulence is 
developed in the normal component, then quantity (5£(uj)5£(—uj)) scales as uj~ 5 / 3 . Interestingly, 
if we accept the conditions of the experiment by Roche et al.[8j], then (e.g. for the mean flow v n ~ 1 
m/s at T = 1.6 K) we have 

(5£(uj)5£(-uj)} w 4 * 10 22 w~ 5 / 3 , 

which is close to the experimental data in order of magnitude. Note also the dependence on the 
applied velocity (about oc v*) is also consistent with experimental data. 

III. HVBK CASE 

In the case of quasi-classical turbulence, the set of vortex line is believed to consist of the 
many bundles containing a large number of threads inside of them. The macroscopic behavior 
of these bundles is quite similar to the dynamics of eddies in ordinary fluids. Th e coarse-grained 
hydrodynamic of the vortex bundles is studied by many authors (see e.g., 



2l|,[2j,|2j,|24|), 
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but the basis of these stidies was the hydrodynamics of rotating superfluids, or the Hall-Vinen 



Bekarevich-Khalatnikov (HVBK) model (see e.g., book 25[]). In the vortex bundles, the vorticity 
field of uj s (to be exact, its absolute value) and the vortex line density £ are related to each other 
with the use of the Feynman's rule: 



(8) 



This formula reflets the fact that the VLD C, in this case coincides with the areal density of lines 
in a plane perpendicular to the bundle. In terms HVBK dynamics of the vorticity obeys to the 
following equation (see [25f |) 



Of 



V x [v L x cj 8 ], 



0) 



where is the velocity of lines, 



vl= v s +a [u) s x (v n - v s )] . (10) 

The meaning of Eq. ([9]) is that it describes the motion of vortex lines in the transverse (with respect 
to the unit vector lu s along the vorticity) direction, when the coarse-grained superfluid velocity v s 
differs from the normal velocity v n . This equation is derived again without fluctuations. To take 
into account the latter, we use the formula ([3|) into the equation (|9|) and have 

dib s (C + 5C) 

— = V x [(v ns + Sv ns )(C + <$£)]. 

After little algebra this equations is transformed to 

dCo s (5£) 



dt 



aCV x 6v ns (11) 



The HVBK theory describes the redistribution of preexisting vortex lines in the transverse direction, 
it does not include a mechanism of the appearance of new lines. In papers [3], [2^], [3] it was shown 
that the bundle structure of quantized vortices develops inside the eddies of the normal component 
due to proliferation of vortex filaments. Mechanism of this proliferation is quite involved, it reminds 
the developing of the vortex tangle in an applied counterflow with the growth of the number lines 
due to reconnections, and with the growth of length due to relative velocity with the subsequent 
polarization. We take into account this process by adding a relaxation-type term in Eq. (fTTTk 
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dCo s {5£) 

at 



a£V x 5v ns H 6£. 




Here t s is the time of relaxation of the bundle structure as described by Samuels et al. [261] . [271] . [2. 
In Fourier the component we have equation 

— ioo8£(u) = ia£—(k x 5v ns ) -\ S£, 

Ps T s 

which lead to the following spectral density: 



Again, as in the case of Vinen turbulence, the shape of the spectrum depends on the value of 
ujt s . The question of the relaxation time is a bit vague and unclear. Analyzing the numerical 
results, Samuels [26] offered an expression for t s , which included the circulation of the vortex 
tube (of normal component) and its size. Therefore it is impossible to apply directly his result. 
The important fact is that quantity r s is proportional 1/y/a, thus, it decreases with increasing 
temperature. 

For large r s , which implies a small coupling between the normal and superfluid components 
and is realized for a small temperature, the spectrum < 5£{uj) 5£(—uj) >oc (5v n (uj)5v n (—uj)) , 
(w -5 ' 3 for Kolmogorov turbulence in the normal component). In this case the result is similar to 
the previous (Vinen turbulence) case, considered above. Even the value of spectrum is close to 
the one, which is induced by the counterflowing turbulence. However, for small ujt s << 1, which 
corresponds to large temperature (strong coupling due to the large mutual friction) the spectrum 
behaves as < 5C{uj) 5£{—uj) >oc uj 2 (<5v n (w)<5v n (— u)), (w 1 / 3 for Kolmogorov turbulence), and the 
intensity of this spectrum is much lower. This result is in good qualitative agreement with the 
numerical result . . ) . As for the quantitative analysis, there is problem of determination 
of the relaxation time r s , whivh is beyond the scope of the present study. 



< 5C(u) 6£(-u) > 



i + K) 2 



(13) 



IV. DISCUSSION AND CONCLUSION 



We studied analytically the spectrum of fluctuations of the vortex line density both in the coun- 
terflowing (Vinen) turbulence and on the basis of the HVBK theory. In both case these deviations 
of quantity 5£{t) appear due to strongly fluctuating field of the normal velocity. Deviations of 5£(t) 
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evolve in the relaxation like manner, which is determined either from the Vinen equation or from the 
HVBK equation. It allows to find Fourier transform and evaluate the spectra < 5£(oj) 5£(—u) >. 
The crucial difference between these two cases is that the stirring force for C(t) in the Vinen case 
is proportional to perturbations of the normal velocity 5v n (t), whereas in the HVBK this force 
is related to V x 5v n (t). This difference results in different spectra and their dependence on the 
temperature. From analysis of the final relations for spectral densities < 8C(cu) 5C{—uj) > we can 
conclude that for the small temperature, both model offer the cj _5//3 specrtrum, whereas for large 
temperature the HVBK theory results in the w 1 / 3 dependence. This conclusion may serve as a 
basis for experimental determination of what kind of turbulence is realized in different types of 
generation. 
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